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ABSTRACT 

O , 

High redshift QSOs provide important opportunities for advancing our understanding of the 
astrophysics of galaxy formation and evolution. Although a growing number of these sources 
are now known to exist at redshifts beyond 4, systematic studies of the emission-line properties 
of these objects are quite limited. We present results of a program of high signal-to-noise spec- 
troscopy for 44 QSOs using the MMT and Keck observatories. The quasar spectra cover 1100 
- 1700 A in the rest frame for sources spanning a luminosity range of ~ 2 orders of magnitude. 
Comparisons between these data and spectra of lower redshift quasars reveal a high degree of 
similarity, although differences are present in the profiles and the strengths of some emission 
features. An examination of the luminosity dependence of the emission lines reveals evidence 
for a weak or absent Baldwin effect among z > 4 QSOs. We compare measurements for objects 
in our sample with results from other high redshift surveys characterized by different selection 
techniques. Distributions of equivalent widths for these different ensembles are consistent with 
a common parent population, suggesting that our sample is not strongly biased, or in any case, 
subject to selection effects that are not significantly different from other surveys, including the 
Sloan Digital Sky Survey. Based on this comparison, we tentatively conclude that the trends 
identified here are representative of high z QSOs. In particular, the data bolster indications of 
supersolar metallicities in these luminous, high-z sources, which support scenarios that assume 
substantial star formation at epochs preceding or concurrent with the QSO phenomena. 

Subject headings: quasars: emission lines — galaxies: abundances — galaxies: evolution 
1. Introduction 
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1 Based in part on observations obtained at the MMT 
Observatory, a joint facility of the University of Arizona The properties of QSOs at the highest known 

and the Smithsonian Institution. redshifts are important for understanding the 

2 Based in part on observations obtained at the W.M. 

Keck Observatory, which is operated as a scientific part- g pace Administration. The Observatory was made possi- 

nership among the California Institute of Technology, the ble by the generous financial support of the W.M. Keck 

University of California and the National Aeronautics and Foundation 
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early history of galaxies and the formation of su- 
permassive black holes. Over the past decade, a 
growing number of quasars have been reported at 
redshifts z > 4, corresponding in look-back time 
to more than ~ 90% of the age of the Universe. 
Objects in this regime are thus inherently young, 
and accordingly of great interest for studying the 
onset of accretion-powered activity in relation to 
galaxy formation and evolution. 

Optical studies of QSOs at z > 4 to date 
have mostly emphasized their discovery, luminos- 
ity function, and use as probes of intervening 
matter seen in absorption. Investigations of the 
emission-line properties of these objects have been 
far more limited (Schneider, Schmidt & Gunn 
1991; Storrie-Lombardi et al. 1996). Further 
scrutiny of their emission spectra in relation to 
those of their low-redshift counterparts is clearly 
desirable. Specifically, such comparisons may re- 
veal information on the timescale of star formation 
and chemical evolution associated with the QSO 
phenomenon (Hamann & Fcrland 1999), and on 
time-dependent behavior of the accretion process 
(Mathur 2000). 

This paper presents a detailed analysis of the 
optical emission line properties for a large sam- 
ple of z > 4 QSOs. High signal-to-noise, mod- 
erate resolution spectra of 44 quasars were ob- 
tained over multiple observing runs at the MMT 
and Keck observatories. The study employs di- 
rect measurements of the emisson lines as well as 
comparisons of composite spectra. Preliminary re- 
sults from this survey were reported by Shields & 
Hamann (1997). Here we compare our spectra 
with those of redshift ~ 2 quasars and find general 
agreement, although subtle differences are present 
in the strengths and profiles of some emission fea- 
tures. Our objects show little evidence for a Bald- 
win Effect, the empirical anticorrclation between 
line equivalent width and luminosity L often found 
in lower redshift samples (Baldwin 1977; see Os- 
mer & Shields 1999 for a review). We also investi- 
gate the degree to which our findings may be influ- 
enced by selection effects. We compare the prop- 
erties of our targets with published measurements 
for other z > 4 QSO samples representing a variety 
of discovery methods. This analysis indicates that 
our sources are representative of known quasars 
in this redshift regime, and suggests that the ob- 
served differences from lower- z sources are at least 



partially intrinsic. In subsequent papers (e.g., Di- 
etrich et al. 2001), we will present more detailed 
comparisons of quasar emission-line properties as 
a function of redshift and luminosity, based on a 
sample of > 800 objects that span a much wider 
range in z and L. 

Although the values of the adopted cosmologi- 
cal parameters play an important role at these ex- 
treme redshifts, we adopt throughout the paper an 
H a = 50 km s -1 Mpc -1 , q a = 0.5, A = cosmol- 
ogy, to ease the comparison with earlier published 
work in this field. 

2. Observations & data analysis 

We began a program of spectroscopy of z > 4 
QSOs in 1994, with targets selected from the 
sources known at that time (~ 50) or reported in 
subsequent years in the literature. Data were ob- 
tained in a series of 7 observing runs at the MMT 
Observatory during 1994 - 1996, and 2 runs at 
the W. M. Keck Observatory in 1997 October and 
1999 June. The goal of this program was to ob- 
tain a statistically significant description of quasar 
spectra in this high redshift regime; 32 objects 
were observed at the MMT, and 12 at Keck, for a 
total of 44 QSOs. Specific targets were selected 
based on their observability with the brightest 
available sources given priority in order to max- 
imize the resulting signal-to-noise ratio. We gen- 
erally avoided sources known to exhibit strong ab- 
sorption lines near the emission redshift. The ma- 
jority of the observed objects are taken from the 
APM Color Survey (Storrie-Lombardi et al. 1996, 
BR/BRI prefix) and the Second Palomar Sky Sur- 
vey (Kennefick et al. 1995b, PSS prefix). Table 1 
lists the redshifts, r magnitudes, and the discovery 
reference for each source. 

Data were acquired with the Red Channel of the 
MMT Spectrograph (Schmidt, Weymann & Foltz 
1989) and with LRIS at the Keck II telescope (Okc 
et al. 1995). Observations were obtained with 
a long slit of 1 "width, yielding full width at half 
maximum (FWHM) spectral resolution of - 10 A 
and ~ 7 A, with pixel sampling of 3.53 A pixel -1 
and 1.86 A pixel -1 , respectively. These resolu- 
tions correspond to 300 - 400 km s _1 across the 
C IV A1549 emission line, or typically 0.07 - 0.1 
times this line's FWHM. The spectral coverage 
was chosen to span the redshifted Lya A1216 - 
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He II A1640 interval; the He II feature is of specific 
interest since the N V A1240/He II A1640 ratio can 
be used as an abundance diagnostic (Hamann & 
Ferland 1993; Ferland et al. 1996). The obser- 
vational setup was essentially identical for the two 
telescopes. A long-pass filter (LP495 at the MMT, 
GG495 at Keck) was used to prevent second-order 
contamination. 

Multiple exposures totaling ~ 1 — 2 hrs were ob- 
tained for each quasar, yielding a typical signal-to- 
noise ratio of ~ 50 per pixel. Observations were 
obtained at an airmass < 1.5, and the slit was 
oriented approximately perpendicular to the hori- 
zon in order to minimize the effects of atmospheric 
differential refraction (Filippcnko 1982). The at- 
mospheric seeing was typically 1" (or better for 
the Keck II telescope), but not all the observa- 
tions were obtained under absolute photometric 
conditions. 

Calibration and reduction of the resulting spec- 
tra were carried out using standard techniques as 
implemented in the IRAF software package 3 . Ob- 
servations of nearly featureless standard stars were 
used to provide flux calibrations and to correct for 
the atmospheric absorption bands. He-Ne-Ar arc- 
lamps were used for the wavelength calibrations. 
Multiple spectra acquired for the same object were 
averaged. Figure 1 presents all the reduced spec- 
tra included in this observational program, sorted 
on RA. 

A small number of known Broad Absorption 
Line (BAL) QSOs were observed as part of the 
sample in order to enable additional compar- 
isons with low-rcdshift quasars. These sources 
are PSS 1048+4407, BR 1144-0723, and PSS 
1438+2538. Our spectra, which have greater S/N 
and wavelength coverage than many of the pub- 
lished discovery spectra, also reveal for the first 
time broad absorption lines in PC 0027+0525 and 
PSS 0137+2837. 

The primary focus of the present work is on 
the emission lines. Accurate measurement of these 
features requires correction for absorption in many 
cases. For the majority of objects, the absorption 
appears as narrow lines due to associated or inter- 

3 The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for 
Research in Astronomy Inc. (AURA), under cooperative 
agreement with the National Science Foundation. 



vening systems, superimposed on either the emis- 
sion lines or continuum. In these cases we removed 
the absorption feature through a simple interpola- 
tion. For the BAL QSOs and several additional 
objects (BRI 0151-0025, PSS 1159+1337, BRI 
1346-0322, PC 1415+3408, and PSS 1435+3057), 
the absorption was severe enough that a reliable 
reconstruction of the unabsorbed spectrum was 
not possible. We excluded these sources from our 
measurements of the emission lines and from our 
composite spectra. All objects in this sample show 
strong Lya forest absorption, which clearly influ- 
ences the measured profile and flux of the Lya 
emission line in many cases; we did not attempt to 
correct for this effect. Objects in our sample are 
subject to small amounts of Galactic foreground 
extinction, spanning Ay = 0.04 — 0.49, with a typ- 
ical value of ~ 0.1 (Schlegcl, Finkbeiner & Davis 
1998). We corrected the spectra for the resulting 
reddening using the empirical selective extinction 
function of Cardelli, Clayton & Mathis (1989). 

The spectra were Doppler-corrected to the rest 
frame using redshift values that we measured for 
the majority of sources from the C IV A1549 line. 
The observed wavelength for each case was estab- 
lished by fitting a Gaussian to the top 20% of the 
profile. For situations where the measured line was 
significantly affected by absorption or low signal- 
to-noise, we verified the result with fits to the top 
50% of the profile. For the BAL QSOs, the red- 
shifts listed in Table 1 are taken from the discov- 
ery reference. For the remaining objects where 
C IV was not measurable due to absorption (see 
above), Lya was used instead. The redshift mea- 
surements in these cases relied on Gaussian fits 
that used as a template only the red side of the 
profile. The resulting z determinations are in good 
agreement with previously published values, with 
an estimated uncertainty of < 0.01. 

Luminosities for our sample targets were de- 
rived using published photometric data. Rest- 
frame 1450A fluxes (AB magnitudes at 1450A) 
were obtained directly from the discovery papers 
for the majority of the objects in our sample. For 
the objects where these numbers were not avail- 
able, the flux measured from our spectra was used, 
after the absolute flux calibration was adjusted 
to obtain consistency with published broad-band 
photometry. The luminosity values, expressed in 
ergs s" 1 Hz -1 , are also included in Table 1. 
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To investigate the emission-line properties of 
the observed QSOs, we constructed composite 
spectra, and employed measurements of the lines 
in individual objects. 

The equivalent widths (EWs) of the strong 
lines, specifically for Lya + N V and C IV A1549, 
were measured interactively, using a direct inte- 
gration of the line flux referenced to the interpo- 
lated continuum level. For the Lya+N V feature 
we used the extension of the continuum from the 
red side only, and did not attempt any deblending. 
The resulting values for C IV are very similar to 
those obtained by modelling the lines with 2 Gaus- 
sians, representing broad and narrow components. 
Possible inaccuracies in these measurements are 
dominated by ambiguities in placement of the con- 
tinuum; experiments with alternative choices of 
the continuum level suggest that this contributes 
a systematic uncertainty of <~ 20%. 

In our analysis, we found it useful to com- 
bine our measurements with published emission- 
line data for other z > 4 QSO samples, and in 
particular those discovered via the Palomar CCD 
Grism Survey (Schneider, Schmidt & Gunn 1991, 
1997; Schmidt, Schneider & Gunn 1995, hereafter 
referred to as SSG, ~ 15 objects), and early dis- 
coveries from the Sloan Digital Sky Survey (SDSS; 
Fan ct al. 1999, 2000, 2001, - 35 objects, only 
those with z > 3.9 were considered in this analy- 
sis). In contrast with our measurements, the pub- 
lished EW values for these samples rely on sin- 
gle component Gaussian fitting. Our results and 
the published EWs can be meaningfully combined 
only if the line measurement procedures are con- 
sistent. We tested different methods of measuring 
the emission lines with our data. For the C IV 
line, the single Gaussian fits give, in general, lower 
EWs than those obtained with our direct integra- 
tion technique, as they tend to lose flux from the 
wings if they are prominent, and from the peak 
if the cores dominate. In order to obtain a rea- 
sonable consistency in the comparisons involving 
the SSG and the SDSS samples, we measured the 
C IV EWs in both ways. For Lya + N V, the 
fit results are sensitive to the method of treating 
the two lines and the asymmetry resulting from 
Lya forest absorption; consequently, we employed 
only the EWs obtained by direct integration of the 
profile. The values are recorded in Table 1. 

Finally, we constructed composite spectra us- 



ing the data for individual sources normalized to 
unit mean flux over the wavelength range 1430 
A- 1470 A. Average and median composites were 
calculated over the wavelength range common for 
all objects. The average composite may be influ- 
enced by several extreme objects with unusually 
strong and narrow lines (BR0401-1711, BRI1050- 
0000, BR2212-1626, BR2248-1242). The median 
spectrum may thus be more representative of the 
typical object in our sample. 

3. Spectroscopic properties 
3.1. General characteristics 

An initial scrutiny of the spectra of z > 4 
QSOs reveals substantial agreement with those of 
their lower redshift counterparts. This similarity 
is noteworthy, given the widely differing amounts 
of time available for the evolution of the AGN 
and underlying host galaxy. In Figure 2 (upper 
panel) we compare our z > 4 average composite 
spectrum with the composite spectrum derived as 
an average 4 from the total Large Bright Quasar 
Survey (LBQS) sample, which is weighted toward 
z w 2 objects in this wavelength range (S. Morris 
1999, private comunication; see Brotherton et al. 
2001 for additional information). This compari- 
son is especially useful given that our sample and 
the subset of the LBQS contributing to this wave- 
length interval have nearly identical average lumi- 
nosities (£i,(1450A) w 31.3), so that any differ- 
ences should reflect a redshift dependence or dif- 
fering selection criteria. We estimated the average 
LBQS L„(1450A) from the absolute B magnitude 
given by Francis et al. (1991), their Figure 8, fol- 
lowing the methodology of Schneider, Schmidt & 
Gunn (1989) to relate the flux level in the two dif- 
ferent bandpasses, while correcting for differences 
in cosmological parameters. In general terms, the 
high and intermediate redshift composites agree 
remarkably well; they have the same strong emis- 
sion lines, similar continuum shape redward of the 
Lya feature, and comparable strengths and pro- 
files of the lines. The pronounced Lya forest in 
the z > 4 composites does not reflect an intrinsic 
difference in the nature of the emission sources, 
but rather an increase in the opacity of the inter- 

4 The algorithm used in constructing the LBQS composite is 
described by Francis ct al. (1991) 
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galactic medium at larger redshifts. 

Similarities between the QSOs at z > 4 and 
those at lower redshift can also be seen in measure- 
ments of spectral features for individual sources. 
Surveys of QSOs at z w 2 - 3 (Francis et al. 1992; 
Francis, Hooper & Impey 1993) display strong 
(negative) correlations of line widths (FWHM) 
with line EWs and with line peak/continuum ra- 
tios. Similar tendencies are present in the high z 
sources. Figure 3 illustrates these results for all 
quasars in our sample, for which the C IV lines 
were measurable. Sources with narrow lines also 
tend to show the strongest emission. Furthermore, 
spectra of QSOs at z < 3 exhibit systematic line 
velocity shifts that correlate with the degree of 
ionization of the emitting species (Gaskell 1982; 
Espey et al. 1989; Tytler & Fan 1992; Mcin- 
tosh et al. 1999). As reported previously by 
Storrie-Lombardi et al. (1996), this behavior is 
also present in z > 4 samples. 

A detailed comparison of the high and interme- 
diate redshift composites in Figure 2 shows excel- 
lent agreement for some lines but dissimilarities 
in others. C II A1335 and Si IV+O IV] A1400 
remain basically unchanged, while differences are 
present in the strengths of Lya, N V, C IV, and 
O I A1304. The Lya emission in z > 4 quasars 
is considerably affected by the forest absorption, 
which may remove up to ~ 50% of the line, but 
the high-z composites nonetheless show a stronger 
Lya feature than in the z m 2 average. The N V 
emission is also enhanced in the high z systems, 
and the same is true for the O I A1304 emission. 
The rest-frame equivalent width of the latter fea- 
ture is ~ 3A in the composite z > 4 spectra, which 
is about 50% larger than that seen in spectroscopic 
surveys at lower redshifts (e. g., Brotherton et al. 
1994; Forster et al. 2000, EW(0 I) « 2A). 

The C IV line is stronger in the z > 4 average 
than in the intermediate z composite spectrum. 
The difference is most pronounced in the line core, 
while the line wings are more nearly comparable 
in the two composites. As noted previously, the 
sample includes several remarkable objects with 
unusually strong and narrow lines, which poten- 
tially distort the comparison when using an aver- 
age composite for this modest sample. To address 
this question, the high z median and the redshift 
~ 2 composite are shown for comparison in Fig- 
ure 2, lower panel. The Lya and C IV features 



are clearly reduced in the high z median, but still 
stronger, especially in the core, than those in the 
LBQS composite. 

Aside from the differences in Lya and C IV, 
the high redshift median and average spectra are 
remarkably similar. In particular, the O I and 
N V emissions are equally strong in both high z 
composites. The degree of variation between the 
members of our sample can be seen in the standard 
deviation of the normalized spectra as a function 
of wavelength, plotted in Figure 4. The largest 
variations are present in the cores of the prominent 
emission lines (Lya and C IV) , consistent with the 
previous results for lower redshift samples (Francis 
et al. 1992; Brotherton et al. 2001). The lower 
panel of Figure 4 shows the standard deviation of 
the mean as a function of wavelength, which pro- 
vides a measure of the uncertainty of the average 
spectrum. 

3.2. O I and N V 

The O I A 1304 and N V A 1240 lines show dif- 
ferences between the intermediate and high red- 
shift spectrum composites that merit further at- 
tention. These features potentially convey infor- 
mation about the structure of the broad-line re- 
gion and the physical conditions in the emitting 
plasma. Additionally, both lines are possible trac- 
ers of metallicity in the QSOs, which is of partic- 
ular interest for understanding these sources in an 
evolutionary framework. 

We first examine whether we have correctly 
identified the O I feature. Previous quasar 
emission line studies have raised the question 
of whether the A 1304 feature is indeed O I 
emission or something else. The O I line is a 
triplet (A1302.17, A1304.86, A1306.03) with a mean 
laboratory wavelength (obtained by weighting 
each transition by its Einstein A coefficient) of 
1303.50A (Verner et al. 1996). For the high op- 
tical depths that may occur in broad line clouds, 
the individual features in the triplet will approach 
equal strengths, modifying the mean to 1304. 36A. 
The feature we attribute to O I is seen at a wave- 
length of ~ 1307A in these composites. At first 
glance, a better identification might seem to be 
with Si II AA1304.37, 1309.27 (mean laboratory 
wavelength 1307. 63A). However, this doublet is 
expected to be accompanied by considerable emis- 
sion in other Si II lines, such as AA1196, 1264, 
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1531, 1817 (Gaskell 1982; Dumont & Mathez 
1981), which are weak or absent. Furthermore, we 
note that the high z composites are derived from 
spectra that are Doppler-corrected using redshifts 
given by C IV, which is blucshiftcd relative to 
the low ionization lines. An identification of the 
A 1307 feature with O I would be consistent with 
this trend, while Si II would not. We conclude 
that the O I A1304 identification is secure. 

The enhanced strength of O I in the high red- 
shift sources may reflect a structural difference in 
their broad- line regions (BLRs). Several studies of 
AGN ensembles suggest that low ionization lines 
(e.g. O I, C II) are emitted preferentially in low- 
velocity BLR components, while Lya and high 
ionization lines (N V, Si IV, C IV) tend to trace 
higher velocity plasma (e. g., Brotherton et al. 
1994; Francis et al. 1992; Wills et al. 1993; Bald- 
win et al. 1996). This is sometimes described in 
terms of a BLR with two components: the Inter- 
mediate Line Region (ILR) and Very Broad Line 
Region (VBLR). In our high redshift sample, the 
narrow component makes a larger fractional con- 
tribution to the spectra than in the lower redshift 
systems. The C IV line has a larger velocity width 
in the LBQS composite (FWHM « 5600 km s" 1 
) than in the z > 4 composites, both average and 
median (FWHM « 4170 km s" 1 and 4900 km s _1 , 
respectively). The unusual strength of O I emis- 
sion may thus be the most obvious tracer of a char- 
acteristically enhanced ILR in z > 4 QSOs. 

The strong N V emission is not expected in sys- 
tems with an enhanced narrow component, how- 
ever, because of its high ionization. The system- 
atically higher N V emission may instead reflect 
high metallicity, consistent with standard chemi- 
cal evolution scenarios for young, massive ellipti- 
cal galaxies or galaxy bulges. Hamann & Ferland 
(1992, 1993) suggested that N V, which is a col- 
lisionally excited line, may be elevated relative to 
other lines as a consequence of the secondary ni- 
trogen enrichment in vigorously star-forming envi- 
ronments. Interestingly, enhancements in the O I 
emission may provide supporting evidence for this 
interpretation. O I A1304 is apparently formed via 
fluorescence with the H Ly/3 line (see, e.g. Netzer 
1990 for a discussion). The strength of the O I 
feature is thus expected to scale with the O/H 
abundance ratio (Netzer & Penston 1976). How- 
ever, the O I line strength may be influenced by 



various other factors, such as the sensitivity of the 
line coincidence to the velocity field in the emit- 
ting region. 

To summarize, the N V and O I features suggest 
that the broad-line regions of QSOs at z > 4 may 
be intrinsically different from those at lower red- 
shift. Low- velocity components of the BLR are 
more prominent in these high-z sources, as sig- 
naled by the strength of O I and the narrow veloc- 
ity width for C IV. In addition to this structural 
difference, however, the strength of N V emission 
implicates characteristically higher metallicities in 
the z > 4 objects. High metallicities may also con- 
tribute to the enhancement of the O I line. 

3.3. The EW vs. Luminosity correlation 

To investigate the amplitude of the Baldwin Ef- 
fect for QSOs at z > 4, we present plots of the rest 
equivalent width of C IV A 1549 as a function of the 
continuum luminosity L v (measured at 1450 A in 
the rest frame), for our observed sample (Figure 5, 
upper panel), and for an enlarged sample that in- 
cludes also measurements obtained from the SSG 
and SDSS surveys (Figure 5, lower panel). The 
dotted line is the relation found by Osmer et al. 
(1994) for a sample of 186 quasars covering a wide 
range of redshift [z < 3.8). The overall results at 
high redshift (z > 4) appear compatible with the 
trend and scatter found for AGNs at lower redshift 
and luminosity, but do not, by themselves, show 
a Baldwin Effect correlation in C IV. A Spear- 
man test is consistent with the null hypothesis at 
> 99% level. The corresponding diagrams for Lya 
+ N V are shown in Figure 6, and likewise are con- 
sistent with no correlation. These measurements 
combine Lya and N V, and therefore, a possible 
trend in the Lya line alone, for which the Baldwin 
Effect has been previously claimed, may be weak- 
ened by the presence of N V, which does not show 
such EW — L v anticorrelation (Osmcr, Porter & 
Green 1994; Korista et al. 1998). 

To explore further the luminosity-dependent 
behavior, we divided the quasars in our sample 
into two luminosity bins, with log L J/ (1450 A) 
= 30.34 - 31.12 and 31.31 - 31.97. Average and 
median composites for these subsets are shown in 
Figure 7. The noise and the intrinsic peculiari- 
ties present in single sources are considerably di- 
minished in the composite spectra, making them a 
valuable tool in identifying line strength variations 



6 



originating in differences in luminosity. A simple 
inspection of the individual lines reveals a small 
decrease in the EW with luminosity for Lya, but 
not for C IV or other emission lines. 

The weak or absent Baldwin Effect in our z > 4 
sample may be an intrinsic characteristic of these 
objects, which are extreme in both redshift and 
luminosity. Alternatively, the absence of a clear 
trend may stem from the limited luminosity range 
of our sample and the substantial scatter that al- 
ways appears in the L v — EW relation. Selec- 
tion effects influencing our sample are also a sig- 
nificant concern; if strong-lined objects are over- 
represented among our high luminosity QSOs, the 
result would be a weakened Baldwin correlation. 
We explore the role of selection effects further in 
the next section. 

3.4. Selection biases 

An important consideration in interpreting the 
properties of our high redshift QSOs is the extent 
to which selection effects can influence the prop- 
erties of a sample. Here we examine the methods 
used for identifying z > 4 QSOs and ways that 
selection criteria may be reflected in emission-line 
characteristics. 

The majority of QSOs known at z > 4 have 
been discovered by color-selection techniques (Ir- 
win et al. 1991; Kennefick et al. 1995a, b; 
Storrie-Lombardi et al. 1996, 2001) that rely in 
particular on the contrast between the observed 
spectral region centered on the Lya feature, and 
the continuum at shorter wavelengths, which is 
strongly absorbed by the Lya forest. Quasars are 
thus distinguished from stars primarily on the ba- 
sis of their extremely red B - R or similar colors, 
which separate them from the stellar locus in the 
color-magnitude/color-color diagrams. The emis- 
sion lines can contribute significantly to the total 
flux observed in typical photometric bandpasses. 
Their role is amplified at these particularly high 
redshifts since the scaling factor (1 + z) boosts the 
values of the observed EWs. 

The prominence of Lya could therefore intro- 
duce biases affecting the detectability of z > 4 
quasars. Most of the objects in our sample come 
from the APM Color Survey (BR/BRI quasars). 
These quasars are drawn from a magnitude- 
limited sample, and as a consequence, the strong- 



lined sources could be over-represented because of 
the contribution of the Lya line to the measured 
flux in the R band (Malmquist bias). The am- 
plitude of this effect can be roughly evaluated by 
estimating the change in the apparent (R) mag- 
nitude produced by the presence of a strong Lya 
line: for a rest frame EWi = 30A (weak line) and 
EW2 = 160 A (strong line) with respect to the 
unabsorbed continuum, the observed (z w 4.5) 
equivalent widths would be EWf s « 165 A and 
EWf « 880 A respectively. With an R band- 
pass of width ~ 1000 A, the line increases the 
measured flux by ~ 16% and 88% compared with 
a pure continuum source. The ratio of the two 
broadband fluxes with strong compared to weak 
Lya would be f2f.f1 — 1-62, which corresponds to 
a difference in magnitude ~ 0.5 mag. This value is 
actually an underestimate of the effect because the 
continuum blueward of Lya is highly suppressed 
by the absorption forest. In addition, for some of 
the observations in the BRI survey, a narrower R 
bandpass (~ 600A) was used, resulting in an in- 
crease of the fractional contribution of Lya to the 
measured flux. This effect leads to a bias in favor 
of objects with large EW, that could potentially 
account for the unusual examples of strong-lined 
QSOs in our sample. The same trend was sug- 
gested previously by Kennefick et al. (1995b), 
who calculated the selection probabilities of QSOs 
(from POSS II observations) for three cases of 
Lya + N V line strength (see their Fig. 6). These 
considerations apply to the BR/BRI, PSS, and 
SDSS sources, which rely on similar bandpasses 
for QSO identification. 

Although this effect seems to be important, it 
could be counteracted by the treatment of the 
other bandpasses in the selection process. If only 
sources that are detected in B are included, as 
is the case for the APM survey, objects with very 
strong lines near the limiting R magnitude may be 
excluded, because their R-band flux is dominated 
by Lya and the continuum (which dominates in B) 
is correspondingly weak. However, the selection 
process is further affected by the fact that, for very 
weak lines, the contrast between the Lya emission 
and the continuum blueward of it may not be high 
enough to satisfy the necessary color criterion. In 
the case of the APM survey, the understanding of 
this bias is complicated by the use of a magnitude- 
dependent color threshold, due to signal-to-noise 
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considerations (Storrie-Lombardi et al. 2001, see 
their Figure 1), again introducing a preference for 
strong-lined sources. The fact that the majority of 
the color selected objects in our sample fall above 
the dotted line in the Baldwin diagram, may be 
due, in part, to this bias. 

A thorough understanding of the Malmquist 
and color biases would require extensive model- 
ing based on the choices of filters, limiting magni- 
tudes, and/or colors. An alternative is to investi- 
gate the EW properties of samples discovered by 
different techniques. We compared the EWs and 
composite profiles of our color-selected sources to a 
grism-selected sample, which includes objects ob- 
served by us and by Schneider, Schmidt & Gunn 
(1991). The SSG sample consists of quasars se- 
lected based on the detection of Lya in slitless 
spectra. As discussed by Schmidt, Schneider & 
Gunn (1995), this detection technique may also 
lead to preferential selection of QSOs with strong 
lines. Figure 8 shows the rest-frame and the ob- 
served EW distributions, for C IV and La, for the 
color-selected and grism-selected samples. The K- 
S (Kolmogorov-Smirnov) statistics are consistent 
with the two samples having the same parent pop- 
ulation. The average and median composite spec- 
tra should indicate similar trends. Figure 9 shows 
this comparison, which reveals a remarkable sim- 
ilarity between the composites. Likewise, there is 
no apparent difference in N V and O I; both lines 
are evidently stronger in the high redshift sources, 
regardless of the selection technique. We can con- 
clude that the color selection methods (primarily 
BR/BRI) find QSOs with analogous properties to 
those of sources found by the SSG grism survey. 
Therefore, the corresponding biases, if apparent, 
affect the samples to a similar degree. 

More useful comparisons may be those with the 
SDSS sample, which is claimed to be less biased. 
According to the authors, the selection probabil- 
ity in the SDSS detection of high redshift ob- 
jects does not strongly depend on the emission 
line strength, but on the redshift and the SED 
shape. Comparisons between observed EW(Lya) 
distributions for our sample, the grism (SSG) sam- 
ple and the SDSS survey show that they are con- 
sistent with a single distribution function (prob- 
ability of being drawn from the same distribu- 
tion is respectively KS sflss/ourS(lm p ie = 0.633, 
KSsDss/grism = 0.910). The rest and observed 



frame EW histograms, for both Lya and C IV, are 
displayed in Figure 8, separately for color-selected 
QSOs from our sample, for the grism-selected ob- 
jects from SSG, and for the SDSS quasars. The KS 
statistic in each case is consistent with a common 
parent population. 

One possible complication in this comparison 
may be that the characteristic luminosities for 
these different samples are not the same. Inspec- 
tion of Figures 5 and 6 indicate that quasars in our 
sample are somewhat more luminous than those in 
the SDSS. The Baldwin Effect would then predict 
smaller EWs for our sources. The fact that the 
EW distributions are actually similar could thus 
be a coincidence, due to a bias favoring strong- 
lined objects in our sample. A closer inspection of 
the luminosity distributions suggests that such a 
conspiracy is unlikely to be significant. There is 
substantial overlap in luminosities for the different 
samples, and differences in their median values are 
small; median L v for our QSOs exceeds that of the 
SDSS sources by only a factor of 2. 

We conclude that the selection effects, if 
present, are not more evident in our sample than 
in the other high redshift surveys. 

4. Discussion 

Interpretations of the strong cosmological evo- 
lution of AGNs are generally based on the func- 
tional form of a (simple) observed luminosity func- 
tion of bright QSOs at high z and prescriptions 
for the growth and luminosity of the underlying 
black holes. In the cold dark matter scenario, 
as a paradigm for a hierarchical cosmogony that 
models the evolution of the luminosity function, 
quasars are believed to be a short-lived, first phase 
of the formation of a galaxy in the potential well 
of a dark matter halo (Hachnelt & Rees 1993; 
Haiman & Locb 2001). 

Observations at high redshifts intercept epochs 
when quasars are very young objects, which thus 
might be expected to have characteristically low 
metallicity. However, the trends reported here 
in the N V and O I broad emission lines (sec- 
tion 3.2) suggest that the metallicities are typically 
higher in higher redshift, and implicitly more lu- 
minous QSOs. The prominent emission lines indi- 
cate that, in spite of their youth, the high-redshift 
QSOs have undergone substantial chemical enrich- 
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mcnt on a timescale that is short compared with 
their life span, implying a rapid star formation 
process in the host galaxy. Therefore, the very ex- 
istence of objects at z > 4 presents severe timing 
problems (time to form and turn on), that can be 
solved only under the assumption that they reside 
in the most massive objects which have collapsed 
at these early epochs, associated with the rare high 
peaks in the primordial density field (Turner 1991; 
Katz et al. 1994). This result may naturally 
arise from the analogy with the mass-metallicity 
relationship present in the low z spirals and el- 
lipticals: massive galaxies reach higher metallici- 
ties because their deeper potential is better able 
to retain their gas against the galactic winds built 
by the thermal pressure from supernovae, while 
low-mass systems eject their gas before high en- 
richments are attained. Thus, quasar metallicities 
should be similarly tied to the gravitational bind- 
ing energy of the local star-forming regions, and 
perhaps also to the total mass of their host galax- 
ies (Hamann & Ferland 1999). 

Consistent with this picture, recent AGN stud- 
ies have revealed correlations between the mass 
of the central black hole and QSO luminosities, 
QSO host masses, and the stellar velocity disper- 
sion (Laor 1998; Magorrian et al. 1998; Ferrarese 
& Merritt 2000; Gebhardt et al. 2000). Quasar 
broad lines trace matter within only the central 
few parsecs, but the emerging black hole/bulge 
connections strongly suggest that the emission-line 
plasma is closely linked to the larger star-forming 
environment. Recent cosmic-structure simulations 
(Gnedin & Ostriker 1997) show that protogalac- 
tic condensations can form stars and reach higher 
than solar metallicities at z > 6, suggesting that 
at epochs preceding, or concurrent with the QSO 
formation period, the bulge star formation already 
occurred. Therefore, interpretation of the abun- 
dance data in these extreme redshift sources yields 
unique constraints on the evolution of those en- 
vironments, indirectly probing the epoch and ex- 
tent of early star formation associated with QSOs. 
The observed enhancements in the N V and O I 
emission reinforce the previous evidence of solar or 
greater metallicities (Z > Z©) that occur before 
QSOs become observable (Dietrich & Wilhelm- 
Erkens 2000; Hamann & Ferland 1999 and refer- 
ences therein). 

The evidence of structural differences reported 



in Section 3.2 for the broad line region in these 
young sources provides further motivation for ex- 
amining their evolutionary context. If the en- 
hanced ILR emission is real and not a consequence 
of selection effects, this behavior may be an impor- 
tant tracer of accretion physics and its relationship 
to redshift and luminosity. Based on our initial re- 
ports of narrow UV lines in z > 4 QSOs, Mathur 
(2000) suggested an analogy between these sources 
and narrow-line Seyfert 1 (NLS1) galaxies. In this 
picture, both classes of objects are in an early evo- 
lutionary phase, in which accretion proceeds at or 
near the Eddington limit. While this scenario has 
several appealing aspects for explaining NLS1 phe- 
nomena, assigning it in general to high-z QSOs 
may be premature. While the latter sources are 
young in cosmological terms (ages < 2 Gyr), this 
does not ensure that they are in an early phase 
of accretion, since the timescale for QSO activ- 
ity in an individual source is estimated to be only 
a few times 10 7 yr (e.g., Kauffmann & Hachnclt 
2000). Furthermore, narrow profiles in the UV 
lines are probably not an adequate basis for link- 
ing our sources to a NLS1 classification, which is 
based on optical lines. As discussed by Wills et 
al. (2000), the UV line widths do not correlate in 
a simple way with H/3 width or with other NLS1 
properties. Further inquiry is needed to under- 
stand the implications of line profile behavior in 
high-redshift QSOs. 

5. Conclusions 

Although an increasing amount of information 
is becoming available on the properties of QSOs 
discovered at very high redshift, only limited ef- 
forts have been made to survey systematically the 
emission- line properties of these objects and/or 
the selection effects related with the techniques by 
which they were discovered. We have conducted a 
thorough analysis of 44 high signal-to-noise spec- 
tra from the MMT and Keck observatories which 
yield measurements of the most prominent emis- 
sion lines in the rest-frame interval llOOA- 1700A. 

Composite spectra for the whole data set and 
for subsets were constructed and analysed in order 
to investigate: a) the emission properties of z > 4 
QSOs in comparison with those of their lower red- 
shift counterparts, b) the luminosity dependence 
of the emission features, and c) the role of selec- 
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tion effects in existing samples. There are several 
conclusions that may be inferred from the work we 
have outlined here: 

1 . In terms of their ultraviolet rest-frame spec- 
tra, z > 4 QSOs strongly resemble quasars at low 
redshift. Subtle differences are present, however, 
and in particular, Lya, N V, C IV and O I are 
stronger in the high z sources. The C IV line is 
also systematically narrower in the z > 4 objects. 

2. Among the high redshift QSOs, a weak Bald- 
win Effect is possibly present in Lya but not in 
C IV or other lines. The lack of a strong trend 
may reflect the limited span of luminosity in the 
existing z > 4 samples. 

3. Selection effects favoring strong-lined objects 
are a significant concern for surveys of high red- 
shift QSOs. Our sample includes several sources 
with extremely strong, narrow and peaked lines 
that may arise from such a bias. However, quanti- 
tative comparisons with other existing surveys of 
z ^ 4 QSOs, including those discovered by SDSS, 
suggest that our sample overall is not strongly bi- 
ased in this way. 

4. All z > 4 composite spectra show strong N V 
and O I. The unusual strengths of these features 
and the narrow C IV profile suggest the presence 
of characteristic structural differences in the BLR, 
and also high abundances of heavy elements (solar 
or up to several times solar metallicities) in quasar 
environments at these early times. These find- 
ings are consistent with standard chemical evolu- 
tion scenarios for young, massive bulge-dominated 
galaxies. 
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Fig. 1. — Observed spectra of QSOs at z > 4. The flux, Fa, is in units of 10 17 erg s 1 cm 2 A : . The 
prominent emission features are marked. 13 
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Fig. 1 continued.- 
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rest-frame wavelength (A) 

Fig. 2. — Composite spectra (average -upper panel, and median -lower panel) of z > 4 QSOs included in this 
observational program are shown in comparison with the LBQS composite, representative of z w 2 objects. 
The y-axis has been normalized to unit mean flux over the wavelength range 1430 A- 1470 A. The redshift 
measured from the C IV feature was used for the Doppler correction of the z > 4 spectra. As a consequence 
of the different velocity shifts presented by the low- and high-ionization lines, the alignment of the high-z 
spectral features with those in the low z composite required a slight additional shift of our composites. 
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Fig. 3. — Equivalent widths and line peak/continuum ratios for C IV emission as a function of C IV line 
widths (full width at half maximum, FWHM), measured in our sample. 



18 



8 



i 1 1 1 1 1 1 1 

standard deviation spectrum 




Fig. 4. — Upper panel: Flux standard deviation in Fa relative to the average composite spectrum, as a 
function of wavelength for the individual spectra comprising the composites. Lower panel: The standard 
deviation of the mean, expressed in %. 
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Fig. 5. — Rest-frame EW of the C IV emission line as a function of absolute luminosity of the quasar (at 
1450 A rest frame). The Baldwin Effect reported by Osmer, Porter & Green (1994), based on measurements 
at lower redshift, is shown by the dotted line. The upper panel shows measurements from our sample only, 
with different symbols for the color (x) and grism (•) selected objects. The lower panel adds measurements 
for z > 4 QSOs from Schneider, Schmidt & Gunn (1991) and the SDSS survey. Although consistent with 
the low redshift trend, the measurements do not, by themselves, show the expected correlation. 
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Fig. 6.— Equivalent width of Lya + N V versus continuum luminosity for z > 4 QSOs. Symbols are the 
same as in Fig. 5. 
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Fig. 7. — Comparison between median and average composites computed for two sets of quasars, grouped 
according to their luminosity at 1450A (L„ expressed in ergs s _1 Hz -1 ). The same normalization as in Figure 
2 is used. A weak drop in EW with luminosity can be noted for Lya but not for C IV or other lines. 
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C IV; three different samples are compared using a KS test criterion which indicates that the distributions 
are consistent with a common parent population. 
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Fig. 9. — Comparison of composite spectra for color-selected and grism-selected QSOs. The composites 
include spectra published by SSG in addition to the data presented here. Average and median spectra are 
shown in the upper and lower panels, respectively. 



24 



Table 1 

qsos observed at mmt and keck telescopes. 



Object 


z a 


r 


b CIV 

EWi c (A) EW 2 d (A) 


Lya + N V logL„(1450) 
EW C (A) 


Rcf 


-TOO UUUO-f-Z/OU 


4. 


,23 


1 Q 

iy 


.3 


33 


32 


95 


31 


,30 f 


Q 
O 


RR (1(11 Q 1 ^99 

i j i v uuiy-iozz 


4. 


.52 


1 Q 

iy 


.0 


43 


34 


97 


31 


.53 


R 



1 V_ UUZ(-(-UOZl 


4, 


,21 


91 

zi 


.9 


57 


43 


133 


30 


,09 


1 3 

10 


1 V. uuz * -t-uozo 


4, 


,09 


91 
Zl 


.5 












1 

10 


jIJOj UUOOZO^rUU 


4. 


,75 


91 
Z 1 


.2 


69 


64 


125 


31 


.11 


A 

4 


roo uuoy+uuuo 


4 


,15 


1 Q 

iy 


.5 


26 


23 


78 


31 


,24 f 


7 
( 


rpt m n^-i-nfYW 


4, 


43 


1 8 


.6 


44 


36 


87 


31 


.52 


1 ^ 


pp ni *?i -urn 
I v_ uioi-t-uiou 


3 


,78 


1 Q 


,4 


17 


15 


72 


31 


.37 


1 n 

1U 


PSJSJ HI '?7-U9S'?7 c 


4, 


,30 


1 Q 

iy 


.0 












Q 



rpt m m nn9 c i 

I J 1 v 1 U101-UUZ0 


4. 


,19 


1 8 


.9 






72 


31 


.44 


1 ^ 
10 


roo UlOZ-t-U/00 


4. 


,05 


1 Q 

iy 


.6 


27 


27 


60 


31 


.06 f 





RPT (19/1 1 MAR 
I J 1 VI UZ41-U140 


4, 


,01 


1 8 
10 


,2 


27 


24 


57 


31 


.62 


1 

10 


PSJSJ (1948-1-18(19 


4. 


.42 


1 X 
lo 


.4 


32 


31 


95 


31 


.76 f 


O 


pp n , ?n7-Ln999 

UoUi-rUZZZ 


4, 


.39 


9(1 
ZU 


,4 


20 


20 


75 


30 


.81 


1 1 
1 1 


RR (14(11 1 71 1 
I J I V U4U1-1 ( ±± 


4, 


.22 


1 8 
10 


,7 


68 


61 


195 


31 


,49 f 


1 R 
10 


rr nam nzL^n 


4, 


.34 


1 s 

10 


.9 


28 


28 


66 


31 


.35 


Q 1 R 

y, io 


RPT (10^9 (11 1 


4, 


.40 


1 R 
lo 


,7 


28 


27 


50 


31 


.55 


1 R 
10 


rrt 1 m c t^nn' : t^ 

I J I V 1 IUIOtuUOO 


4, 


.38 


1 & 
10 


,8 


36 


35 


56 


31 


.43 


1 A 
14 


rr 1 n c t97 

I J IV 1UOO-UOZ ( 


4. 


.49 


1 8 
lo. 


.5 


26 


24 


72 


31 


.53 


R '\R 
O, 10 


pqq 1 n/is_L/i/in7 c 


4, 


.45 


1 Q 

iy 


.3 












a 
o 


rrt 1 n^n nnnn 

1 J I V 1 ±uou-uuuu 


4, 


.28 


1 R 
10 


,7 


73 


73 


119 


31 


.26 


1 A 
14 


roo 1U0 1 -{-4000 


4. 


.10 


1 7 
1 ( 


,7 


24 


21 


74 


31 


,79 


a 
o 


RRT 1 1 (IS (17ZL7 


3. 


.91 


1 8 
10 


.1 


28 


25 


93 


31 


.44 


1 R 
10 


rrt 1 1 1 n^ni of; 

I J I V 1 111U I U1UO 


3. 


.91 


1 8 

10 


.3 


20 


19 


42 


31 


.40 


1 R 
IO 


RR 1 1 zLzL n79 c ! c 


4. 


16 


1 R 
io 


.6 












1 R 
io 


-Too noy^ioo/ 


4. 


.08 


1 s 

10 


.5 






15 


31 


,91 f 


■3 



RR 1909-0795 


4. 


.69 


1 8 


.7 


17 


17 


22 


31 


.87 


fi 1 6 

1U 


PC 1301+4747 


3. 


.99 


21 


.3 


46 


46 


27 


30 


.32 


12 


PSS 1317+3531 


4. 


.36 


19 


.1 


30 


29 


85 


31 


,39 f 


8 


BRI 1328-0433 


4. 


.20 


19. 


.0 


39 


26 


71 


31 


.37 


9, 16 


BRI 1346-0322 


4. 


.00 


18. 


.8 






104 


31 


.21 


16 


PC 1415+3408 


4, 


.58 


21 


.4 






70 


30 


.66 


13 


PSS 1435+3057 


4. 


.35 


19 


.3 






20 






7 


PSS 1438+2538 6 


4. 


.24 


19 


.5 












8 


PC 1450+3404 


4. 


.19 


20 


.8 


44 


38 


102 


30 


.68 


13 


BRI 1500+0824 


3 


.93 


18. 


.8 


23 


18 


44 


31 


.32 


15 


PSS 1618+4125 


4. 


.21 


19 


.6 


31 


26 


73 


31 


,06 f 


3 


Q 1745+6226 


3. 


,89 


18 


.8 


38 


34 


79 


31 


,80 f 


1 


RX 1759+6638 


4. 


,32 


20 


.0 


45 


40 


68 


31 


.30 


5 


Q 2203+2900 


4. 


,41 


20 


.8 


31 


26 


53 


30 


.87 


2 



25 



Table 1 — Continued 



Object 


z a 


r b 
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IV 




Lya + N V 


logL„(1450) 
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11 



Calculated from our spectra using the C IV emission-line 

b photometry in r or analogous bandpass (R for the APM survey, r for PSS, r* for SDSS, r± for 
SSG) as reported in the respective discovery papers 

c rest-frame EWs calculated by direct integration of the line flux 

d rest-frame EWs calculated by fitting the line to a single Gaussian 

c exhibits broad absorption lines (BAL); z values are taken from the discovery papers 

f values derived using the flux fx (1450) measured from our spectra; estimated errors are approx- 
imately ±0.2 dex. 

References.— (1) Becker et al. (1992); (2) Dickinson & McCarthy 1987; (3) S. Djorgovski, 
private comunication; (4) Fan et al. 1999; (5) Henry et al. 1994; (6) Isaak et al. 1994; (7) 
Kennefick et al. 1995a; (8) Kennefick et al. 1995b; (9) Kcnncfick et al. 1996; (10) Schmidt et al. 
1987; (11) Schneider et al. 1989; (12) Schneider et al. 1991; (13) Schneider et al. 1997; (14) Smith 
et al. 1994a; (15) Smith et al. 1994b; (16) Storrie-Lombardi et al. 1996. 
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